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FOREWORD 


The  theory  of  thermally  activated  radiationless  transitions, 
which  compete  with  luminescence,  is  well  understood  in  principle, 
but  there  is  a  paucity  of  examples  of  its  detailed  quantitative 
application  to  well  characterized  systems,  especially  in  the 
intermediate-coupling  regime  occupied  by  transition-metal 
complexes.  An  analogous  situation  prevails  in  the  case  of 
excited-state  absorption.  Athermal  fluorescence-quenching 
mechanisms,  which  operate  in  the  strong-coupling  regime  occupied 
by  color  centers,  are  more  controversial.  The  aggressive 
extension  of  theory  has  been  made  urgently  imperative  by  the 
exploitation  of  such  systems  in  potential  tunable-laser 
materials,  where  fluorescence  quenching  and  excited-state 
absorption  are  major  loss  mechanisms.  Low-crystal -field  chromium 
complexes  in  ordered  perovskites  of  cubic  elpasolite  structure 
are  a  paradigm  for  vibronic  laser  materials  based  on  transition- 
metal  doping  of  ionic  crystals.  The  laser-active  T1®(1)  color 
center  in  KCl  and  its  analogues  exemplify  the  strong-coupling 
regime.  The  reported  investigation  has  helped  to  identify 
factors  controlling  these  important  loss  mechanisms  by  detailed 
experimental  characterization  of  the  systems  in  question  and  by 
successful  application  of  semi-empirical  models.  Ab-initio 
computer  modeling  by  embedded-molecular-cluster  calculations, 
which  require  further  refinement,  may  ultimately  facilitate 
formulation  of  design  criteria  for  such  materials. 
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STATEMENT  OF  PROBLEM  STUDIED 


Fluorescence  quenching  and  excited-state  absorption  are  major 
loss  mechanisms  in  potential  tunable-laser  materials.  The  object 
of  the  reported  investigation  was  to  achieve  a  quantitative 
understanding  of  the  factors  which  control  such  processes.  The 
ultimate  goal  is  to  develop  theoretical  models  with  genuine 
predictive  capability,  in  order  to  provide  guidance  in  the 
selection  of  potential  tunable-laser  materials  by  permitting 
their  a-priori  evaluation  with  respect  to  these  loss  mechanisms. 
The  subject  of  the  present  report  on  ARO  Contract  No. 
DAAL03-89-K-0059  is  the  completion  of  an  investigation  which  was 
initiated  under  ARO  Contract  No.  DAAG29-82-K-0158  and  continued 
under  ARO  Contract  No.  DAAL03-86-K-0017 .  Experimental  techniques 
employed  in  this  investigation  for  characterizing  low-field 
chromium  compl&^es  and  for  evaluating  radiationless-transition 

J 

rates  include  photoluminescence  and  Raman  spectroscopy,  at  both 
ambient  and  high  pressures,  and  two-photon  excitation 
spectroscopy.  Radiationless-transition  rates  were  calculated  by 
computer  line-shape  simulations.  A  semi -empirical  model 
constrained  by  spectral  information  served  to  identify  the 
dominant  photoluminescence  quenching  mechanism.  Ab-initio 
embedded-molecular-cluster  calculations  and  lattice  dynamics 
calculations  were  employed  in  the  development  of  theoretical 
models  with  potential  predictive  capability. 
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SUMMARY  OP  RESULTS 


I.  Introduction 

Low  crystal -field  complexes  of  Cr^*  are  characterized  by  broad¬ 
band  (^T2g  -»  ^A2g)  fluorescence,  in  contrast  with  high-field 
complexes  which  exhibit  narrow-band  (^Eg  -*  ^A2g)  phosphorescence. 
This  property,  together  with  a  large  Stokes  shift,  make  these 
low-field  complexes  ideal  candidates  for  tunable  solid-state 
vibronic  laser  appl  ications .  However,  the  strong 
electron-lattice  coupling  which  lends  useful  properties  to  these 
and  analogous  materials  also  enhances  their  susceptibility  to 
both  thermal  quenching  of  fluorescence  and  excited-state 
absorption,  related  fundamental  processes  which  proceed  from  the 
relaxed  excited  state  of  the  laser-active  center. 

The  project  reported  here  involved  a  coordinated  experimental  and 
theoretical  investigation  of  thermal  fluorescence  quenching  and 
excited-state  absorption  in  low-crystal-field  chromium  complexes 
in  ordered  perovskites  of  cubic  elpasolite  structure.  This  class 
of  materials  possesses  the  virtue  that  a  trivalent  cation 
impurity  can  be  accomnodated  in  a  rigorously  octahedral  site 
without  charge  compensation.  Additionally,  it  provides  chemical 
sequences  with  a  range  of  lattice  parameters  and  crystal  fields. 
The  specific  materials  chosen  for  study  include  K2NaGaP(:Cr^^, 
K2NaScFg :  Cr^*  and  CS2NaYCl5:Cr^*. 
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Laser-active  color  centers  exhibit  even  stronger  electron-lattice 
coupling  than  transition-metal  complexes,  making  them  susceptible 
to  an  athermal  fluorescence  quenching  mechanism  which  proceeds 
from  the  unrelaxed  excited  state.  The  final  phase  of  this 
investigation  included  a  theoretical  investigation  of  the  Tl**(l) 
center  in  KCl  and  its  analogues  which  is  still  in  progress. 

Under  the  initial  contract.  No.  DAAG29-82-K-0158 ,  fluorescence 
spectra  and  lifetimes  of  chromium- doped  elpasolites  were  measured 
as  functions  of  temperature  in  order  to  parameterize  models  and 
to  determine  activation  energies  and  pre-exponential  frequency 
factors  for  radiationl ess-transition  rates. ^  A  semi-empirical 
theoretical  model  was  proposed  in  which  radiationless 
deactivation  of  the  chromium  complex  is  mediated  by  quadratic 
coupling  to  an  asymmetric  vibrational  mode.^  Photoluminescence 
spectra  of  the  fluoride  elpasolites  measured  in  a  diamond-anvil 
cell  revealed  a  pressure-induced  re-ordering  of  excited  states, 
with  a  consequent  transition  from  low-field  to  high-field 
behavior,  and  a  pressure-induced  enhancement  of  frequencies 
associated  with  defect-related  vibrational  modes. 

Both  experiments  and  computer  modeling  were  extended  during  the 
period  of  the  subsequent  contract.  No.  DAAL03-86-K-0017 . 
Photoluminescence  lifetimes  in  the  chloride  elpasolite,  measured 
as  functions  of  temperature  and  pressure  in  the  diamond-anvil 
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cell,  revealed  a  pressure-induced  shift  of  the  onset  of  thermal 
luminescence  quenching  to  higher  temperatures.^  Ab-initio 
embedded  cluster  calculations  on  the  ground  electronic  state  of 
the  chromium  complex  yielded  predictions  of  the  pressure 
dependence  of  local  compressibility  and  defect-related  vibration 
frequencies  in  good  accord  with  experiment.  Raman  scattering 
measurements  in  the  diamond-anvil  cell  provided  information  about 
the  pressure  dependence  of  both  defect-related  and  host-lattice 
vibration  frequencies . 

Substantial  additional  progress  in  both  theory  and  experiment  was 
made  during  the  period  of  the  contract  covered  by  this  report. 
No.  DAAL03-89-K-0059 .  Photoluminescence  spectroscopy  was 
extended  to  higher  pressures  and  temperatures  with  the 
construction  of  a  new  diamond-anvil  cell,  revealing  an  excited- 
state  level  crossing  in  the  chloride  elpasolite  at  higher 
pressure  than  in  the  fluorides,  and  additional  Ixuninescence 
lifetime  measurements  on  the  chloride  elpasolite  revealed  that 
the  pressure  dependence  of  the  thermal  activation  energy  for 
fluorescence  quenching  is  in  excellent  agreement  with  predictions 
of  the  semi-empirical  model. Extension  of  embedded- cl  us  ter 
calculations  to  the  excited  state  provided  optical  transition 
energies  in  excellent  agreement  with  experiment Radiationless 
transition  rates  predicted  by  optical  line-shape  simulations 
derived  from  embedded-cluster  calculations  agreed  less  well  with 
experiment .  Lattice-dynamics  calculations  provided  a  good 
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account  of  Raman  spectra^^'^^  and  of  line  shapes  associated  with 
resolved  vibronic  structure  of  emission  spectra.  Two-photon 
excitation  spectroscopy  provided  information  about  electron 
lattice  coupling  in  a  higher  excited  state, Theoretical  models 
were  developed  to  explain  anomalous  features  of  two-photon 
spectra.  Computer  modeling  of  laser-active  color  centers  by  ab- 
initio  embedded  cluster  calculations,  initiated  during  the  period 
of  this  contract,  and  is  still  in  progress;  preliminary  results 
are  encouraging.  These  developments  are  elaborated  and  assessed 
in  the  following  sections. 

II.  Thermal  Luminescence  Quenching  Mechanism 

There  exists  an  extensive  literature  on  the  theory  of  thermally 
activated  radiationless  transitions  spanning  more  than  fifty 
years.  In  the  Mott  theory, the  activation  energy  for  thermal 
fluorescence  quenching  is  identified  with  the  difference  between 
the  crossing  energy  of  the  adiabatic  potential  energy  surfaces 
and  the  minimxim  energy  of  the  surface  associated  with  the  excited 
electronic  state.  The  modern,  formal  development  of  the  theory 
proceeds  from  the  pioneering  contributions  of  Huang  and  Rhys^^  and 
of  Lax^^.  Since  no  radiationless  transitions  are  possible  between 
exact  eigenstates  of  a  physical  system,  the  radiationless 
transition  rate  depends  on  the  non-stationary  state  which  is 
prepared  in  a  given  experiment,  and  that  in  turn  depends  on  the 
coherence  time  of  the  exciting  radiation. In  the  adiabatic 
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coupling  scheme,  radiationless  transitions  are  assumed  to  proceed 
from  Born-Oppenheimer  states  and  are  mediated  by  the  non- 
adiabaticity  operator.  First-order,  time-dependent  perturbation 
theory  then  leads  to  an  expression  for  the  radiationless 
transition  rate  Wj^  of  the  form^^ 

WjTO  =  (AJp[nG(Qo+Up)  +  (n+l)G(Qo-(Op)  ]  ,  (1) 

n  =  [exp(hUp/kBT)  -  1]'^.  (2) 

It  is  convenient  to  distinguish  between  promoting  modes,  which 
mix  electronic  states  of  different  symmetries,  and  accepting 
modes,  which  absorb  their  energy  dif f erence .  Coupling  to 

accepting  modes  is  reflected  in  the  normalized  line-shape 
function  G(Q),  where  Q  is  measured  from  the  zero-phonon  line  and 
hOo  is  the  energy  gap  to  be  spanned  in  the  multi-phonon 
transition.  The  promoting-mode  interaction  is  incorporated  in 
the  factor  t*Up,  where  (jp  is  the  promoting-mode  angular  frequency. 

Several  investigators  have  demonstrated  the  equivalence  of 
adiabatic-  and  static-coupling  schemes  at  a  certain  level  of 
approximation, and  have  discredited  the  popular  "Condon 
approximation"  introduced  by  Huang  and  Rhys.  The  static-coupling 
scheme  provides  the  most  tractable  formula^®  for 

=  (n/hwp2)  i<#(£;0)iaH/aQpi^(i;0)>|2,  O) 
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where  Qp  is  the  promot in^-node  configuration  coordinate. 

The  central  problem  in  the  application  of  radiationless- 
transition  theory  is  evaluation  of  the  normalized  line-shape 
function  G(Q) ,  in  Eq .  (1),  which  also  describes  the  wavelength 

dependence  of  the  fluorescence  spertrum  for  values  of  its 
argument  0  which  are  very  far  from  the  range  of  values  accessible 
to  direct  observation.  The  line-shape  function  is  given  by 

G(Q)  =  ^  Pa  ^  []  l<Xfp3lX,a,>l^  6(Ofp-Oi„+Oo-Q)  ,  (4) 

a  ^  a 

Pa  =  exp(-hQia/'‘BT)  /  ^  exp( -hOj y/lcgT)  ,  (5) 

1 

where  <XfPjlXia,^  an  overlap  integral  of  vibrational  wave 

functions  for  accepting  mode  a,  associated  with  final  (f)  and 
initial  (i)  electronic  states,  with  quantum  numbers  ^  and  a 
respect i vel y . 

Photoluminescence  lifetimes  of  low-field  chromium  complexes, 
plotted  as  functions  of  temperature,  display  a  characteristic, 
nearly  discontinuous  change  in  slope  which  signals  the  onset  of 
thermal  quenching.  Analysis  of  lifetime  data  for  the  three  host 
crystals  studied  in  the  present  investigation,  KjNaGaFgtCr^^, 
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KjNaScFgtCr^'*’  and  CsjNaYClgtCr^'*',  yielded  thermal -quenching 
activation  energies  of  9240,  7270  and  4250  cm~^ ,  respectively, 
and  pre-exponential  frequency  factors  of  order  10*^  sec"^.  The 
measured  temperature  dependence  of  photoluminescence  spectra 
served  to  characterize  the  chromium  complex. 

The  challenge  during  the  initial  phase  of  this  project  was  to 
account  for  the  observed  thermal -quenching  behavior  in  terms  of  a 
semi -empi r i cal  theoretical  model  whose  parameters  were  severely 
constrained  by  spectral  information.  The  simple  Mott  model 
described  above  failed  by  twenty  orders  of  magnitude  to  explain 
the  observed  radiationless  transition  rate. 

In  application  of  the  more  formal  theory  encapsulated  in  Eqs . 
(l)-(5),  a  promoting  mode  of  tj^  symmetry  is  required  in  order  to 
couple  the  ^T2g  excited  state  with  the  ground  state. 
Counter-rotating  displacements  of  anion  octahedra  about  trivalent 
and  monovalent  cations  were  found  to  provide  a  strongly  coupled 
mode  of  the  required  symmetry, contrary  to  earlier  speculation 
concerning  promot ing-mode  selection  rules. The  promot i ng-mode 
factor  Uf  calculated  from  Eq.  (3)  with  free-ion  chromium  d 
orbitals  and  the  exact  nuclear  potential  for  the  ligands,  is  of 
the  order  of  10*^  sec'*  for  all  three  host  crystals. 

It  was  found  that  required  values  of  the  accept ing-mode  line- 
shape  function,  G(Qo±(i)p),  are  extremely  sensitive  to  the  precise 
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physical  model  adopted  for  the  system  when  the  parameters  of  the 
model  are  constrained  by  empirical  spectral  information  and  by 
theoretical  estimates  of  the  promot i n^-mode  factor  u.  Models 
based  on  linear  couplings  to  modes  of  aj^  symmetry  were  found  to 
fail  by  six  to  eight  orders  of  magnitude  to  explain  the  observed 
r ad i at i on  1  ess- 1 r ans i t i on  rate,  even  when  a  range  of  vibration 
frequencies  was  considered. 

The  failure  of  linear-coupling  models  suggested  consideration  of 
quadratic  coupling;  i.e.,  different  vibration  frequencies  in  the 
two  electronic  states.  Although  quadratic  coupling  to  ai^  modes 
is  possible  in  principle,  it  is  ruled  out  in  the  present  system 
by  direct  experimental  evidence.  The  octahedral  chromium  complex 
has  additional  degrees  of  freedom,  however,  and  coupling  to  lower 
symmetry  accepting  modes  is  permitted  by  the  orbital  degeneracy 
of  the  ^Tjg  excited  state.  In  order  to  limit  the  number  of 
adjustable  parameters,  a  simplified  model  was  considered  which 
includes  just  two  electronic  states,  the  ground  state  and 
one  orbital  component  of  the  ^Tjg  excited  state,  with  linear 
coupling  to  a  single  ajg  accepting  mode  and  quadratic  coupling  to 
a  single  tjg  accepting  mode.  The  single  a^g  accepting  mode  in 
this  model  is  an  effective  mode  representing  both  a|g  and  eg 
modes,  which  actually  have  comparable  linear  coupling. 
Similarly,  the  single  tjg  accepting  mode  in  this  model 
incorporates  all  of  the  quadratic  coupling  which  is  actually 
distributed  between  eg  and  tjg  modes.  Normalized  line-shape 
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functions  were  calculated  for  each  node  by  the  aethod  of  Struck 
and  Fonger,^^  in  which  individual  vibrational  overlap  integrals 
are  evaluated  by  the  Manneback  recursion  fornula,^^  and  the 
composite  line  shape  was  calculated  by  numerical  convolution. 
Both  the  ajg  and  ground-state  t2g  vibration  frequencies  were 
treated  as  adjustable  parameters^  but  the  strength  of  quadratic 
coupling  was  calculated  ab  initio  from  a  point-ion  model.  Not 
only  was  it  possible  to  fit  the  temperature  dependence  of 
fluorescence  spectra  and  the  activation  energy  for  thermal 
quenching  -simultaneously,  but  the  model  also  yielded  the  correct 
order  of  magnitude  for  the  absolute  radiationless-transition  rate 
when  the  estimated  value  of  the  promoting  mode  factor  was  used. 
Furthermore,  the  model  worked  equally  well  for  all  three 
elpasolite  compounds  considered,  and  the  adjusted  ground-state 
tjg  vibration  frequencies  are  in  good  agreement  with 
independently  determined  values. 

In  the  most  recent  phase  of  this  project,  the  semi -empi r i cal 
model  was  applied  successfully  to  the  pressure  dependence  of 
thermal  fluorescence  quenching.  A  model  for  thermal  fluorescence 
quenching  based  on  ab-initio  embedded  cluster  calculations  was 
also  introduced.  These  more  recent  developments  are  considered 
in  subsequent  sections  of  this  report. 

III.  Pressure  Dependence  of  Chromium  Photoluminescence 
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Photoluminescence  experiments  under  intense  hydrostatic  pressure 
were  initiated  to  assist  in  the  development  and  testingT  ot 
physical  models.  In  the  previous  phases  of  this  project, 
pho  t  o  I  umi  nescence  spectra  and  lifetimes  of  K2NaGaFg :  Cr^^  and 
KjNaScFg : Cr^^  were  recorded  with  chopped  He-Cd  laser  excitation 
at  pressures  between  0  and  6.5  GPa  and  at  controlled  temperatures 
between  150  and  300  K,  by  means  of  a  diamond-anvil  cell  mounted 
in  a  cold-finger  dewar.  In  both  compounds  a  transition  from  low- 
field  to  high-field  behavior,  reflecting  a  crossing  of  the 
and  ^Eg  excited  energy  levels,  is  evidenced  by  a  pronounced  blue 
shift  of  the  featureless  fluorescence  spectrum  and  its  ultimate 
replacement  by  a  highly  structured  phosphorescence  spectrum.  The 
interpretation  of  the  vibronic  structure  of  the  high-pressure, 
1 ow- t emperatur e  phosphorescence  spectra  was  accomplished  by 
analogy  with  the  ambi ent -pressure ,  1 ow-temperature  spectrum  of 

the  high-field  compound  KjNaAlFgtCr^^,  and  confirmed  by  Raman 
scattering.  Vibration  frequencies  were  found  to  increase  by  ai 
2%/GPa  in  both  compounds.  The  vibronic  structure  in  the  scandium 
compound  exhibited  an  additional  feature:  the  zero-phonon  line 
and  each  replica  were  split,  providing  evidence  for  a  low- 
temperature  phase  transition  which  was  subsequently  confirmed  by 
both  Raman  and  electron  spin  resonance  spectroscopy. 

Photoluminescence  spectra  of  CsjNaYClgtCr^^  were  measured  with 
Ar-ion-laser  excitation  at  pressures  between  0  and  11.5  GPa  and 
at  controlled  temperatures  between  80  and  300  K.  This  compound 
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exhibits  an  even  more  pronounced  blue  shift  than  the  fluoride 
compounds  and  its  excited-state  level  crossing  occurs  at  a  much 
higher  pressure  (8  GPa,  compared  with  1  GPa  for  the  gallium 
compound  and  3  GPa  for  the  scandium  compound).  There  is  evidence 
for  a  pressure- induced  phase  transition  above  9.5  GPa. 

The  observed  blue  shift  and  level  crossing  in  all  three  compounds 
is  attributed  to  a  substantial  pr essur e- i nduced  enhancement  of 
the  crys t al -f i el d  parameter  Dq,  associated  with  reduction  of  the 
lattice  parameter.  A  pressure-dependent  local  compressibility 
extracted  from  the  measured  blue  shift  of  Cs2NaYCl5:Cr^^  is 
plotted  in  Fig.  1. 

Luminescence  lifetimes  of  Cs^NaYClgtCr^^  measured  as  functions  of 
temperature  and  pressure  are  plotted  in  Fig.  2.  Excitation  was 
provided  by  the  second  harmonic  of  a  Q-switched  Nd-YAG  laser  and 
s i gnal -averaged  data  were  recorded  by  a  LeCroy  fast-transient 
digitizer.  Lifetime  data  were  fit  to  an  expression  of  the  form 

T‘l  =  To"*  exp(oT)  +  s  exp[-(AEo  +  nP)/kBTlt  (6) 

where  the  first  term  on  the  right-hand  side  represents  the 
radiative  transition  rate  and  the  second  term,  the  non-radiat  i  ve 
rate.  The  parameter  values  tq"*  =  6370  sec"*,  a  =  0.00403  K"* ,  s 
=  3.8  X  10*5  sec'*  and  AEo  =  4250  cm"*  were  established  at 
ambient  pressure.  The  activation  energy  for  thermal  fluorescence 
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Fig.  1.  Local  coapressibi  1  i  ty  in  CsjNaYClc ;Cr^^:  (a) 
calculated  for  ground  state  (solid  line)  and  first  excited  state 
(broken  line)  and  (b)  froa  experiaental  blue-shift  data 


Fig.  2.  Photoluainescence  lifetiae  t  of  Cs2NaYCl(:Cr^^  vs. 
teaperature  at  0  GPa  (  ^ )»  0.4  GPa  (  O  )  *nd  1.0  GPa  (  ^  ) 
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quenching  was  found  to  depend  linearly  on  pressure  with 
coefficient  n  =  1668  ±  52  cia'VGPa. 

The  semi -empi r i cal ,  t wo-accept i ng-mode  model  described  in  the 
preceding  section,  which  incorporates  simultaneous  linear 
coupling  to  an  a]^  mode  and  quadratic  coupling  to  a  t^g  mode,  had 
been  devised  to  explain  ambient  pressure  lifetime  data.  With  the 
pressure  dependence  of  its  parameters  derived  from  that  of  the 
photoluminescence  spectrum,  this  model  predicts  t}  -  1758  cm'VcPa 
for  the  pressure  dependence  of  the  activation  energy  in 
Cs2NaYClg:Cr^''’,  in  excellent  agreement  with  experiment. 

High-pressure  phot o 1 umi nescence  measurements  on  chromium 
complexes  in  elpasolites  have  proved  to  be  particularly  revealing 
of  their  properties  and  processes.  Both  the  pressure- induced  re¬ 
ordering  of  excited  energy  levels  and  the  pressure  dependence  of 
thermal  quenching  are  dramatic  results  which  challenge 
theoretical  interpretation.  In  particular,  the  latter  result 
provides  convincing  confirmation  of  the  semi -emp i r i cal ,  two- 
accepting  mode  model  for  thermal  fluorescence  quenching. 

IV.  Coherent  Raman-Induced  Radiationless  De-Excitation 

At  the  suggestion  of  Dr.  Mikael  Ciftan,  an  experiment  to  verify 
by  direct  observation  the  participation  of  specific  accepting 
modes  in  the  quenching  of  chromium  fluorescence  was  attempted  in 
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collaboration  with  Professor  Walter  Bron  at  the  Irvine  canpus  of 
the  University  of  California.  The  technique  of  tine-resolved 
coherent  anti-Stokes  Raman  scattering  (TRCARS)^^  was  adapted  to 
that  purpose.  The  object  of  the  experiment  was  to  accelerate 
radiationless  de-excitation  of  substitutional  chromium  ions  in 
the  cubic  elpasolite  KjNaScFj :  Cr^'*’  by  pulsed  excitation  of 
accept i ng-mode  vibrations.  Dual  synchronously  pumped  and 
synchronously  amplified  mode-locked  dye  lasers,  with  their 
frequency  difference  tuned  to  the  vibration  frequency  of  the  aj^ 
defect  mode,  effected  temporally  coherent  excitation  of  that 
mode.  The  dye  lasers  also  served  to  excite  the  chromium  ions  to 
their  ^Tjg  electronic  state.  The  intention  was  to  monitor 
radiationless  de-excitation  by  detecting  anti-Stokes  Raman 
scattering  of  a  delayed  probe  beam  by  the  accept ing-mode 
vibrations  coherently  excited  in  the  multi-phonon  decay  process. 
The  relative  amplitudes  of  the  several  accepting  modes  would  then 
indicate  the  extent  of  their  participation. 

Phase  matching  was  irrelevant  in  this  experiment  since  the 
impurities  lacked  translational  symmetry;  however,  the  consequent 
loss  of  spatial  coherence  raised  the  threshold  for  the  desired 
effect  by  many  orders  of  magnitude.  The  minimum  energy  per  pulse 
required  to  stimulate  radiationless  de-excitation  was  estimated 
to  be  2.5  joules/cm^,  approximately  coincident  with  the  laser- 
damage  threshold,  as  are  most  other  non-linear  optical  effects 
involving  point  defects.  Unfortunately,  when  the  experiment  was 
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attempted,  laser  damage  occurred  before  any  evidence  of  the 
desired  effect  could  be  detected.  It  was  regretfully  concluded 
that  the  method  may  be  feasible,  but  not  in  this  material. 

V.  Computer  Modeling  of  Low-Field  Chromium  Complexes 

Computer  modeling  of  pressure-dependent  optical  properties  by  ab- 
initio  embedded-cluster  calculations  was  initiated  with  the 
ultimate  objective  of  developing  models  with  genuine  predictive 
capability.  Mo  1  ecu  1 ar -o r bi t a  I  calculations  were  performed  on  a 
twenty-one  atom  cluster,  AgB0CrXg*^^,  for  each  compound 
AjBMXgtCr^^,  by  means  of  an  RHF-LCAO-SCF  program  called  MELD.^® 
The  Is,  2s  and  2p  shells  of  chromium  and  its  nearest  chlorine 
neighbors  were  replaced  by  effective  core  potentials,  while 
double-C  quality  basis  sets  for  valence  orbitals  were  explicitly 
included.  All  nearest-neighbor  fluorine  orbitals  were  included, 
but  the  nearest  sodium  and  potassium  ions  were  represented  by 
bare  effective  core  potentials. 

Interactions  of  ions  inside  the  cluster  with  those  outside  were 
approximated  by  central,  rigid-ion  pair  potentials  of  Buckingham 
form. 


V(r)  =  ZjZjeVr  +  Aj jexp(-r/Pi j)  -  Cjj/r®.  (7) 

Simultaneous  relaxation  of  the  molecular  cluster  and  its 
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surrounding  lattice  was  accomplished  by  means  of  a  modification 
of  the  HADESIII  lattice  statics  code,^^  called  HADESR,  developed 
under  the  present  project.  Electronic  structure  calculations 
were  performed  for  both  the  ^Ajg  ground  state  and  the 
excited  state  for  each  of  the  compounds  of  interest.  The  total 
energy  of  the  cluster  was  calculated  in  each  state  on  a  5x5x5 
mesh  of  symmetrical  displacements  of  X~,  A'^  and  ions 

surrounding  the  chromium  impurity.  This  function  was  fit  to  a 
polynomial  in  symmetry-adapted  coordinates  and  was  used  to  update 
the  total  energy  and  forces  at  each  iteration  of  the  lattice- 
statics  calculation.  The  calculation  then  yielded  a  fully 
symmetrical  relaxed  configuration  of  the  crystal,  including  the 
embedded  cluster,  for  each  electronic  state. 

Asymmetric  displacements  were  treated  more  approximately  by 
minimizing  the  total  energy  of  the  ^T^^  state  with  respect  to  eg 
and  tjg  displacements  of  X~  ions  alone,  with  the  remaining  ions 
fixed  in  their  symmetrically  relaxed  configuration.  The 

correction  for  asymmetrical  relaxation  of  the  remaining  ions  was 
explored  at  ambient  pressure  and  found  to  be  small. 

The  total  energy  of  each  state  was  calculated  in  the  relaxed 
configuration  of  the  other  state  in  order  to  determine  optical 
transition  energies  in  accordance  with  the  Frank-Condon 
principle.  The  effects  of  pressure  were  simulated  by  varying  the 
lattice  parameter,  which  serves  as  a  boundary  condition  by 
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constraining  the  ionic  configuration  at  large  distances  froa  the 
impur i ty . 

The  predicted  -»  optical  emission  energy  of 

K2NaScFg :  Cr^'*'  is  plotted  as  a  function  of  pressure  in  Fig.  3. 
Also  shown  for  comparison  are  the  peak  wavelengths  of  the  broad 
optical  emission  bands  measured  in  a  diamond-anvil  cell.  The 
recorded  emission  spectra  were  corrected  not  only  for  the 
response  of  the  detection  system,  but  also  for  the  frequency 

dependence  of  the  transition  probability,  in  order  to  facilitate 
comparison  with  theory.^®  It  is  evident  that  the  predicted 

emission  energies  are  in  excellent  agreement  with  experiment. 

The  predicted  ^Ajg  -♦  ^Tjg  optical  absorption  energy  of 
K2NaScFg :  Cr^*  at  ambient  pressure,  14,840  cm’*,  also  compares 

favorably  with  the  experimental  value,  15,601  cm'*.  However,  the 
difference  between  the  absorption  and  emission  energies,  the 

Stokes  shift,  is  more  sensitive.  The  calculated  value,  1,771 
cm'*,  increases  to  2,045  cm'*  when  corrected  for  Jahn-Teller 
coupling  to  the  t2g  mode  and  for  eg  relaxation  of  the  remaining 
ions  outside  the  octahedron,  but  it  still  falls  somewhat  short  of 
the  experimental  value,  2,748  cm'*.  The  corrected  Stokes  shift 
predicted  for  Cs2NaYCl5:Cr^^,  1,941  cm'*,  also  falls  short  of  the 
room-temperature  experimental  value,  2,642  cm'*,  but  actually 

exceet-j  a  reported  low-temperature  value,  1,740  cm'*. 
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KjNaScFerCf 


Fiff.  3.  Pressure  dependence  of  the  photol u«i nescence  peak 
energy  for  KjNaScFj !Cr^*.  Theory:  - ,  Experiment:  A. 


Fig.  4.  Pressure-dependent  ground-state  vibration 

frequencies  for  K2NaGaF|:Cr**.  Theory:  - 1  (a)  ai|t  (*>)  Cf 

(c)  t}..  Experiment:  O  (au)t  A  <e|)  and  x  (tjg) 


Approximate  vibration  frequencies  associated  with  the  chromium 
impurity  were  calculated  as  functions  of  pressure  by  exploring 
the  dependence  of  total  energy  on  aj^,  e^  and  tjg  symmetry- 
adapted  displacements  of  halogen  ions  with  respect  to  the  relaxed 
configuration  in  each  electronic  state.  The  calculated  pressure 
dependence  of  ground-state  vibration  frequencies  of  KjNaGaFg ; Cr^'*' 
is  compared  with  values  inferred  from  resolved  vibronic  structure 
of  emission  spectra  in  Fig.  4.  Calculated  vibration  frequencies 
exceed  experimental  values  by  &(  15  %  in  all  three  compounds;  this 
systematic  discrepancy  is  attributed  in  part  to  neglect  of 
electron  correlation. 

A  new  method  was  developed  for  optical  line-shape  simulation  by 
direct  d i agona 1 i za t i on  of  the  matrix  of  the  vibrational 
Hamiltonian  within  each  electronic  state  in  a  common  harmonic 
oscillator  basis.  The  vibrational  overlap  integrals  in  Eq.  (4) 
were  evaluated  as  the  direct  products  of  eigenvectors,  and  the 
corresponding  energy  eigenvalues  were  employed  in  evaluating  the 
delta-function  factors.  This  method  accurately  replicates  line 
shapes  for  combined  linear  and  quadratic  coupling  which  were 
calculated  previously  by  the  Manneback  recursion  relations. 
However,  the  new  method  is  much  more  general  in  that  it  can  also 
accommodate  any  form  of  anharmoni c i ty  for  which  the  matrix 
elements  can  be  calculated;  i.e.,  any  which  is  expressible  as  a 
power  series  in  the  configuration  coordinate. 
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The  new  method  was  applied  to  calculation  of  radiationless 
transition  rates  between  ^Tjg  and  states  of  the  chromium- 

doped  elpasolites  KjNaGaFg,  KjNaScFg  and  CsjNaYClg.  Adiabatic 
potential  energy  surfaces  for  each  electronic  state  were  derived 
from  the  ab-initio  embedded  cluster  calculations  described 
previously.  The  5x5x5  mesh  provided  not  only  linear  and 
quadratic  coupling  to  a^g  modes,  but  also  cubic  and  quartic 
enharmonic  coupling.  Only  linear  and  quadratic  coupling  to  eg 
and  tjg  modes  was  considered,  except  for  the  chloride  compound, 
where  anharmonic  coupling  to  a  tjg  mode  was  investigated  by  means 
of  conf  i  gur  a  t  i  on- i  nt  er  ac  t  i  on  (Cl)  v.a  t  cu  1  a  t  i  ons  . 

Radiationless  tra..sition  lates  predicted  for  KjNaScFg :  Cr^'*'  are 
compared  with  e>  leriment  in  Fig.  5.  Rates  predicted  for  the 
fluoride  compounds  in  the  harmonic  a^>^roximat  i  on  were  lower  by 
twr  orders  magnitude  than  the  values  inferred  from  experiment; 
<^quival  r  <  the  temperature  predicted  for  the  onset  of 
luminescence  quenching  was  too  high  by  a  factor  of  two.  This 
result  is  a  substantial  improvement  on  the  single-configuration- 
coordinate,  linear-coupling  model  with  parameters  constrained  by 
empirical  spectral  information,  which  fails  by  seven  orders  of 
magnitude.  The  predicted  transition  rate  for  the  chloride 
compound,  on  the  other  hand,  was  too  small  by  ten  orders  of 
magnitude,  essentially  because  the  anticipated  quadratic  coupling 
to  the  tjg  mode  failed  to  materialize  in  that  case.  Inclusion  of 
ajg  anharmoni c i ty  worsened  the  agreement  in  all  three  compounds. 
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Fi*.  5.  Radiationless  rate  Wki(T)  for  KiNaScFerCr^*  at  P  = 
0;  (a)  Experiaent*  (b)  Haraonic  theory  and  (c)  Anharaonic  theory 


Fig.  6.  Calculated  (  ■  ■■■■—)  and  aeaaured  Raaan  frequencies  of 
CsjNaYClc  vs.  pressure:  (a)  a|,,  (b)  e^,  (c)  tj})  and  (d)  tjj^ 
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VI.  Lattice  Dynanics  Calculations 


Host-lattice  vibration  frequencies  were  calculated  as  functions 
of  wave  vector  from  the  same  phenomenological  pair  potentials 
employed  in  lattice-statics  calculations.  This  investigation  led 
to  a  re-evaluation  of  the  pair  potentials  subsequently  employed 
in  embedded  cluster  calculations.  Revised  pair  potentials 
incorporate  theoretically  determined  Van  der  Waals  coefficients 
and  a  slightly  reduced  (<  5  %)  Born-Mayer  parameter  p  adjusted  to 
fit  solid-state  properties  at  ambient  pressure. 

The  frequencies  of  Raman-active  modes  at  zone-center  were 
calculated  as  functions  of  hydrostatic  pressure.  Calculated 
frequencies  for  Cs2NaYClg  are  compared  with  experimental  values 
inferred  from  Raman  spectra  in  Fig.  6. 

Vibration  frequencies  were  calculated  along  symmetry  directions 
in  the  first  Brillouin  zone  at  ambient  pressure,  yielding  phonon 
dispersion  curves.  An  imaginary  frequency  associated  with  a  t]^ 
displacement  explains  the  low-temperature  phase  transition 
observed  in  KjNaScFg, 

Frequencies  were  also  calculated  as  functions  of  pressure  on  a 
Kellerman  mesh  in  the  first  Brillouin  zone,  yielding  the  phonon 
density  of  states  and  frequency-dependent  per f ect- 1  at t i ce  Green's 
function  matrix  elements  for  symmetry-adapted  displacements  of 
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the  six  fluorine  neighbors  of  a  trivalent  cation.  The  latter 
were  employed  in  conjunction  with  force-constant  changes  derived 
from  embedded  cluster  calculations  on  the  chromium  impurity  to 
construct  defect-lattice  Green's  function  matrix  elements.  The 
predicted  density  of  normal  modes  projected  on  halogen 
displacements  is  consistent  with  the  narrow  vibrational  sidebands 
in  optical  spectra  which  are  conventionally  attributed  to  normal 
modes  of  the  chromium  hexahalide  octahedron. 

VII.  Two-Photon  Excitation  Spectroscopy 

Two-photon  excitation  spectroscopy  at  low  temperature,  employed 
previously  in  transition-metal  complexes  by  McClure  and  co- 
workers,^*^  has  the  virtue  that  two-photon  transitions  between 
electronic  states  of  the  same  parity  are  electric-dipole  allowed, 
whereas  one-photon  transitions  are  forbidden.  Consequently,  the 
vibronic  structure  of  the  two-photon  spectrum  is  greatly 
simplified,  since  progressions  of  even-mode  phonon  replicas  are 
built  on  a  prominent  zero-phonon  line  rather  than  on  false 
origins  associated  with  odd  modes.  The  object  of  the 
investigation  reported  here  was  to  characterize  the  electron- 
lattice  coupling  in  a  higher  excited  state  of  K2NaScFg;Cr^'*’  in 
order  to  achieve  a  better  understanding  of  excited-state 
absorption  (ESA).  In  particular,  it  is  important  to  account  for 
the  exceptional  width  of  ESA  bands  in  transition-metal  complexes, 
which  makes  them  so  deleterious  to  laser  operation.  Experiments 
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were  performed  on  an  oriented  single  crystal  at  low  temperatures 
(4-10  K)  with  Raman-shifted,  Nd : YAG-pumped  dye-laser  excitation. 
Two-photon  absorption  (TPA)  was  monitored  by  photomultiplier 
detection  of  polarized  chromium  fluorescence. 


The  ^Ajg  -»  ^Tjg,  TPA  transition  was  emphasized  in  the  present 
investigation.  The  corresponding  ESA  band,  ^Tjg  -♦  ^Tjg,, 

overlaps  only  a  small  portion  of  the  tuning  range  of  a  potential 
tunable  laser  based  on  this  material.  However,  the  more 
important  ^Tigj,  state  is  inaccessible  by  the  present  technique 
because  its  TPA  spectrum  overlaps  the  one-photon  ^Ajg  -»  ^T2g 
excitation  spectrum  for  chromium  fluorescence. 


A  we  1 1  - r eso 1 ved  two-component  zero-phonon  line  was  observed  whose 
splitting  was  attributed  to  the  established  1 ow- t emper a t ur e  phase 
transition  in  this  material;  spin-orbit  splitting  is  evidently 
suppressed  by  strong  Jahn-Teller  coupling.  Convolved  line-shape 
simulations  established  equal  coupling  to  a]g  and  eg  modes  in  the 
^T]g,  state,  comparable  in  strength  to  that  for  the  ^Tjg  state, 
plus  much  stronger  coupling  to  the  tjg  mode,  which  is  virtually 
absent  in  the  ^Tjg  state.  The  latter  coupling,  together  with 
opposing  displacements  of  the  eg  coordinate  in  the  two  states, 
accounts  very  well  for  the  width  of  the  excited-state  absorption 
spectrum. 

Additional  features  encountered  in  the  ^Ajg  -»  ^T]ga  TPA  spectrum 
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include  anomalous  polarization  anisotropy  of  the  broad  vibronic 
sideband,  relative  to  that  of  the  zero-phonon  line,  and 
anomalously  extended  progrress  i  ons  of  phonon  replicas.  The  former 
feature  was  explained  in  terms  of  departures  from  the  crude  Born- 
Oppenheimer  approximation  induced  by  strong  coupling  to  the  t^^ 
mode.  The  latter  feature,  which  had  been  observed  previously  in 
another  material,  was  explained  by  the  coupling  of  intermediate 
states  with  diffuse  wavef unct ions  to  lattice  modes  involving 
displacements  of  more  remote  ions. 

The  -*  ^Tjg  TPA  spectrum  was  also  observed,  in  spite  of  the 
fact  that  this  transition  is  forbidden  by  two-photon  selection 
rules. No  zero-phonon  line  was  observed,  and  it  finally  became 
apparent  that  the  entire  spectrum  consisted  of  progressions  built 
on  a  false  eg  origin.  Lineshape  simulations  supported  equal 
coupling  to  a]g  and  eg  modes,  in  confirmation  of  one-phonon 
spectra,  and  negligable  coupling  to  the  tjg  mode.  Neither 
anomalous  polarization  anisotropy  nor  anomalously  extended 
progressions  were  observed.  However,  the  broad  sideband  was 
perturbed  by  Fano  ant i -resonances  associated  with  ^Eg  and  ^T]g 
multiplicity-forbidden  TPA  transitions. 

VIII.  Computer  Modeling  of  Laser-Active  Color  Centers 

Vibronic  lasers  based  on  color  centers  in  ionic  crystals  are 
emerging  as  serious  candidates  for  long-wavelength  sources  in 
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telecommunications.  However,  their  very  strong  e 1 ec t ron- I  at t  i  ce 
coupling  makes  them  susceptible  to  an  additional,  athermal 
fluorescence  quenching  mechanism.  This  quenching  mechanism, 
which  operates  during  the  relaxation  process  at  T  =  0,  was  first 
identified  by  Dexter,  Klick  and  Russell  (DKR)^^  and  was 
subsequently  elaborated  by  Bartram  and  Stoneham*^  and  others.  The 
DKR  rule  predicts  radiationless  de-excitation  of  color  centers 
when  the  excited-state  energy  reached  in  a  Frank-Condon 
transition  from  the  ground  state  lies  above  the  intersection  of 


ad i aba  t i c 

potential  energy 

curves . 

For  1  inear 

coupl ing , 

this 

criterion 

corresponds  to  A 

>  0.25, 

where  A  is 

def i ned  as 

the 

ratio  of 

the  excited-state 

lattice 

relaxation 

energy  to 

the 

optical  absorption  energy. 

The  T1<>(1)  center  in  alkali  halides,  produced  by  subjecting  a 
thallium-doped  crystal  to  ionizing  radiation  and  consisting  of  a 
substitutional  neutral  thallium  atom  adjacent  to  an  anion 
vacancy, has  proved  to  be  one  of  the  most  promising  laser-active 
color  centers. The  electronic  and  vibronic  structure  of  this 
center  were  investigated  in  some  detail  by  multiple  resonance 
t echn i ques . The  analogous  Ga^(l)  and  In^(l)  centers  in  KCl , 
however,  exhibited  little  or  no  luminescence  even  at  a 
temperature  of  10  K.  It  was  demonstrated  by  a  model  calculation 
that  the  disparate  behavior  of  these  analogous  centers  could  be 
explained  in  terms  of  the  DKR  quenching  mechanism,  and  that  the 
parameter  which  distinguishes  them  is  their  spin-orbit  coupling. 
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This  analysis  led  to  the  prediction  that  the  analogous  Pb''’(l) 
center  in  halides  with  divalent  cations  should  exhibit  efficient 
luminescence,  at  least  by  the  DKR  criterion.  That  prediction  has 
since  been  corroborated  for  Pb-doped,  X-irradiated  alkaline-earth 
fluorides, where  laser  action  is  inhibited  by  excited-state 
absorpt ion,^^  and  lasing  was  demonstrated  in  X-irradiated,  Pb- 
doped  KMgF3  which  is  probably  attributable  to  this  center. 

Ab-initio  embedded  cluster  calculations  on  Tl^(l)  centers  were 
initiated  as  part  of  the  present  project.  A  double-C  quality 
basis  set  with  effective  core  potentials  for  the  thallium  atom 
and  its  nearest  chlorine  neighbor  opposite  the  vacancy  was 
augmented  by  diffuse  vacancy-centered  orbitals  derived  from  an 
independent  F-center  calculation.  The  remaining  ions  in  the  44- 
ion  cluster  were  represented  by  bare  effective  core  potentials. 
RHF-SCF-LCAO  molecular-orbital  calculations  were  performed  on  the 
ground  state  and  lowest  excited  state  by  means  of  the  MELD 
program. 

The  spin-orbit  splitting  t  was  reduced  from  that  of  the  free 
thallium  atom  by  the  fraction  of  thallium  p  orbital  in  the  open- 
shell  molecular  orbital.  The  energies  of  the  three  lowest  states 
were  then  approximated  as  follows: 

=  -“2“  *  ^  ]  (A  -  l)*  +  87*  +  Eq.  (8) 
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Wj  =  A  +  7  +  Eq, 


(9) 


E-  +  2Ee 

Eo  =  - - - ,  (10) 


where  E^  and  E^  are,  respectively,  the  energies  of  the  non¬ 
degenerate  ground  state  and  the  doubly-degenerate  lowest  excited 
state,  obtained  from  molecular-orbital  calculations  in  the 
absence  of  spin-orbit  interaction.  The  parameter  yq  was 
introduced  to  correct  for  the  difference  between  the  axial 
chlorine  neighbor  with  valence  orbitals  and  the  transverse 
chlorine  neighbors  represented  by  effective  core  potentials. 


Molecular-orbital  calculations  were  performed  for  various 
combinations  of  axial  thallium  displacements  and  four  symmetry- 
adapted  displacements  of  the  immediate  neighbors  of  the  thallium 
atom  and  chlorine  vacancy.  Simultaneous  relaxation  of  the 
embedded  cluster  and  surrounding  lattice  was  accomplished  by  the 
HADESR  program.  Optical  transition  energies  and  the  ground-state 
thallium  vibration  frequency  are  compared  with  exper iment^^'^^  in 
Table  1. 
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Table  1.  Predicted  and  Measured  optical  transition  energies 
and  thalliuB  vibration  frequency  of  the  Tl^(l)  center  in  KCl 


E,bs  (eV) 

Ee,  (eV) 

u  (10^^  sec"*) 

Theory 

1.38 

0.98 

0.77 

Exper iment 

1.20 

0.83 

0.57 

Adiabatic  potential  energy  functions  generated  in  this  type  of 
calculation  will  ultimately  provide  the  basis  for  a  theoretical 
investigation  of  the  athernal  luminescence  quenching  behavior  of 
inactive  Ga^(l)  and  In*^(l)  analogues,  and  thus  for  a  prediction 
of  their  residual  quantum  efficiency  of  fluorescence. 

IX.  Summary  and  Conclusions 

Radiationless  transitions  generally  belong  to  one  of  two 
categories:  static  processes  which  are  thermally  activated  from  a 
metastable  state,  and  dynamic  processes  which  occur  during  rapid 
relaxation  iaunediately  following  excitation.  Both  categories 
were  of  interest  in  the  context  of  the  present  investigation  of 
luminescence  quenching,  the  former  for  chromium  complexes  and  the 
latter  for  color  centers. 

The  chromium-doped  elpasolite  crystals  selected  for  this 
investigation  provided  a  chemical  sequence  of  low-field  chromium 
complexes,  with  rigorous  octahedral  site  syauaetry,  typical  of 
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transition-metal  based  tunable  laser  materials.  The  range  of 
crys t al -f i el d  parameters  was  further  extended  by  application  of 
hydrostatic  pressure,  culminating  in  the  crossing  of  excited- 
state  energy  levels.  Detailed  characterization  of  these 
materials  was  accomplished  by  one-  and  two-photon  optical 
absorption  spectroscopy,  pressure-dependent  photoluminescence  and 
Raman  spectroscopy  and  electron  paramagnetic  resonance 
spectroscopy.  Photoluminescence  lifetime  data  as  functions  of 
temperature  and  pressure  provided  a  measure  of  thermal 
fluorescence  quenching,  and  the  temperature  and  pressure 
dependence  of  pho t o 1 umi nescence  spectra  served  to  constrain  the 
parameters  of  various  e 1 ec t ron- 1  at t i ce  coupling  models. 

A  semi -empi  r  i  c<t !  ,  two-accept  ing-mode  model,  devised  to  explain 
photoluminescence  lifetime  data  at  ambient  pressure,  proved  to  be 
extremely  successful  in  predicting  the  observed  pressure 
depe.idence  of  the  activation  energy  for  thermal  fluorescence 
quenching  as  well.  Although  the  semi -empi r i cal  model  is  known  to 
be  oversimplified,  this  successful  prediction  appears  to  confirm 
its  essential  features:  the  dependence  of  the  quenching  mechanism 
on  quadratic  coupling  and  on  coupling  to  asymmetric  modes. 

An  ab-initio  three-accept ing-mode  model,  based  on  RHF-LCAO-SCF 
embedded  cluster  calculations  with  static  lattice  relaxation, 
predicts  pressure-dependent  optical  transition  energies  in 
excellent  agreement  with  experiment  in  all  three  compounds.  The 
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predicted  ambient-pressure  Stokes  shifts  are  systematically  30  % 
lower  than  experimental  values  measured  at  room  temperature; 
however,  that  for  the  chloride  compound  is  10  %  higher  than  a 
reported  low-temperature  value.  Predicted  pressure-dependent 
def ec t -r el  a t ed  vibration  frequencies  are  systematically  high  by 
about  15  %.  Finally,  predicted  radiationless  transition  rates 
are  low  by  two  orders  of  magnitude  in  the  fluoride  compounds,  or, 
equivalently,  the  predicted  temperature  for  the  onset  of  thermal 
quenching  is  too  high  by  a  factor  of  two.  This  result  is  a 
substantial  improvement  on  the  s i ng 1 e- conf i gur a t i on-coord i na t e , 

1 i near -coup  1 i ng  model  with  parameters  constrained  by  empirical 
spectral  information,  which  fails  by  seven  orders  of  magnitude. 
However,  the  predicted  rate  for  the  chloride  compound  is  too  low 
by  ten  orders  of  magnitude.  Anharmoni c i ty  worsened  the  agreement 
in  all  cases . 

Both  the  semi -empi r i cal  model  and  the  ab-initio  model  developed 
in  this  project  go  far  beyond  what  has  been  attempted  previously 
in  the  theory  of  radiationless  transitions.  Models  incorporating 
linear  and/or  quadratic  coupling  to  a  single,  syauaetrical  mode 
have  been  reported,  but  they  have  invariably  incorporated 
adjustable  parameters  whose  adjusted  values  are  grossly 
incompatible  with  the  optical  properties  of  the  system  in 
question.  A  proposed  explanation  of  high  radiationless 
transition  rates  based  on  anharmonic  coupling, which  is  widely 
accepted,  has  been  shown  to  be  inapplicable  in  the  present  case. 
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The  semi -emp  i  r  i  ca  1  model  served  to  identify  the  controlling: 
mechanism  for  radiationless  de-activation  of  the  chromium 
complexes  under  investigation.  The  ab-initio  model,  on  the  other 
hand,  still  falls  short  of  the  kind  of  predictive  capability 
required  for  tunable-laser  design  applications,  although  it 
represents  a  giant  step  in  that  direction. 

An  obvious  refinement  of  the  ab-initio  embedded-cluster  model 
would  be  the  addition  of  configuration  interaction  (Cl)  at  the 
level  of  single  and  double  excitations  selected  by  diagrammatic 
perturbation  theory.  The  principal  difference  in  treatment  of 
the  chloride  compound  and  the  fluoride  compounds  was  the  use  of 
an  effective  core  potential  for  the  chloride  ion  and  an  all¬ 
electron  calculation  for  the  fluoride  ion;  the  all-electron 
representation  of  the  chloride  ion  is  feasible  as  well,  but 
pushes  comput ing-t ime  limits.  Even  with  these  refinements, 
however,  the  model  as  presently  formulated  would  be  incapable  of 
predicting  the  observed  temperature  dependence  of  the  Stokes 
shift,  which  is  related  to  both  thermal  expansion  and  quadratic 
coupling.  The  model  portrays  the  lattice  as  too  rigid  at  the 
elevated  temperatures  where  radiationless  transitions  actually 
occur,  and  thus  underestimates  the  el ect ron- 1  at t i ce  coupling  on 
which  the  transition  rate  depends  very  sensitively.  In 

principle,  lattice  relaxation  at  a  finite  temperature  should 
minimize  the  total  Gibbs  free  energy  rather  than  the  total 
internal  energy.  It  remains  to  incorporate  such  a  refinement  in 
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the  basic  lattice-statics  alg^or  i  thin. 


Lattice  dynamics  calculations  served  to  demonstrate  that  the  same 
pair  potentials  employed  in  the  lattice-statics  calculations 
could  explain  pressure-dependent  Raman  scattering  data  as  well. 
They  also  confirmed  the  assignment  of  resolved  vibronic  structure 
of  emission  spectra  to  modes  of  the  chromium  hexahal ide 
octahedron,  and  thus  validated  the  more  approximate  method  of 
calculating  vibration  frequencies. 

Two-photon  excitation  spectroscopy  revealed  strong  Jahn-Teller 
coupling  to  tjg  modes  in  the  state  which  is  absent  in  the 
^Tjg  state,  thus  elucidating  the  origin  of  the  exceptional  width 
of  the  ^T2g  -♦  excited-state  absorption  band.  Unexpected 
and  theoretically  challenging  features  of  the  TPA  spectra 
included  anomalous  polarization  anisotropy,  anomalously  extended 
progressions  and  forbidden  transitions. 

Preliminary  embedded-cluster  calculations  with  static  lattice 
relaxation  performed  on  low-symmetry  laser-active  color  centers 
successfully  predicted  optical  transition  energies  and  a  Raman 
frequency.  These  calculations  are  ultimately  directed  toward  an 
in-depth  investigation  of  dynamic  radiationless  transitions  which 
impede  laser  action  in  analogous  centers. 
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